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•^he  primary  focus  of  research  during  this  funding  period  has  been  to 
obtain  information  with  the  highest  signal-to-noise  ratio  from  faint 
astronomical  sources  using  a  charge-coupled  device  (CCD).  The  starting  point 
of  the  investigation  was  to  determine  the  optimal  sampling  function  for  CCDs 
considering  the  noise  spectrum  characteristic  of  real  devices;  the  noise 
spectrum  of  most  CCDs  can  be  accurately  described  by  a  combination  of  1/f  and 
white  noise.  Then,  a  CCD  system  was  constructed  so  that  all  operating 
parameters  could  be  easily  adjusted.  Some  unexpected  improvements  in 
performance  were  obtained  by  optimizing  the  operating  temperature  and 
increasing  the  number  of  electrode  voltages  which  could  be  independently  set 
compared  to  other  CCD  systems.  ✓The  next  phase  of  the  investigation  was  to 
develop  an  image  processing  system  to  extract,  in  a  statistically  correct  way, 
as  much  information  as  possible  from  astronomical  images.  The  system  was  then 
applied  to  problems  ranging  from  speckle  interferometry  to  global  cosmological 
tests  including  an  investigation  to  determine  the  nature  of  the  dark  matter  in 
the  halos  of  spiral  galaxies. 


The  Final  Report 


Research  during  the  funding  period  of  this  grant  can  be  divided  into  three 
main  areas:  the  development  of  a  charge-coupled  device  (CCD)  imaging  system, 
the  development  of  software  to  analyze  CCD  astronomical  images,  and  the 
scientific  progress  obtained  both  with  and  without  the  CCD.  We  shall  begin  by 
providing  a  brief  introduction  to  CCDs.  The  emphasis  in  our  discussion  will  be 
to  provide  information  which  might  be  of  interest  to  anyone  building  a  CCD 
system  and  which  may  not  be  generally  known. 

A  CCD  is  a  mosaic  of  photodetectors  fabricated  on  a  silicon  substrate. 
CCDs  represent  an  enormous  step  forward  in  imaging  detectors  in  that  they  have 
properties  which  are  close  to  that  of  an  ideal  detector.  CCDs  have  quantum 
efficiencies  of  50  per  cent,  are  linear  over  dynamic  ranges  of  10,000  and  have 
stabilities  of  better  than  one  part  in  1000.  While  imaging  systems  based  on 
CCDs  appear  to  have  enormous  potential,  it  has  been  very  difficult  to  build 
working  systems  which  exploit  the  device's  potential.  The  electronics  to 
operate  CCDs  are  complex  and  there  are  many  operating  parameters  which  must  be 
carefully  adjusted  for  the  device  to  function,  much  less  to  function  well. 

One  of  our  primary  accomplishments  has  been  to  work  out  several  basic 
issues  related  to  CCD  signal  extraction.  We  have  calculated  the  optimal 
sampling  function  for  CCDs  so  that  information  can  be  extracted  as  effectively 


as  possible.  Also,  we  have  determined  the  signal-to-noise  ratio  for  the 
optimal  sampling  function  so  that  it  is  possible  to  compare  several  popular 
sampling  functions  with  the  optimal  to  know  what  improvements  might  be 
possible.  These  results  have  been  derived  for  CCDs  which  exhibit  a  combination 
of  1/f  and  white  noise.  Such  noise  spectra  closely  model  the  noise  in  real 
CCDs. 


Using  these  theoretical  results  as  a  guide  we  built  a  very  flexible  CCD 
system  so  that  virtually  all  CCD  parameters  could  be  easily  and  conveniently 
adjusted  to  optimize  the  response  of  the  system.  To  accomplish  this,  all  29 
timing  parameters  can  be  changed  from  a  computer  keyboard  in  increments  of  100 
ns. 
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In  addition,  we  found  it  was  important  to  keep  the  electronics  as  flexible 
as  possible.  Our  system,  designed  for  a  thinned  RCA  CCD,  was  built  with  17 
adjustable  power  supplies.  There  is  some  advantage  in  using  this  number  of 
supplies  over,  for  example,  11,  but  the  advantage  is  not  large.  Using  less 
than  11  supplies,  however,  would  result  in  degraded  performance.  Specifically, 
we  found  that  it  was  very  important  to  have  separate  high  and  low  voltage 
controls  for  each  of  the  three  horizontal  clocks.  Also,  we  found  that  the  CCD 

performance  was  improved  by  using  three  different  voltage  levels  for  the 

vertical  clocks.  On  the  other  hand,  we  found  that  separate  supplies  for  each 
of  the  three  vertical  clocks  were  not  particularly  important.  Other  supplies 
are  for  phase  reset  and  for  RD. 

During  the  process  of  optimizing  the  performance  of  the  CCD  we  were 
surprised  to  find  that  the  resolution  of  the  CCD  depended  sensitively  on  the 
operating  temperature  of  the  device.  There  was  an  optimal  temperature  with  a 
characteristic  width  of  about  10  degrees  Centigrade.  We  have  no  way  of  knowing 
whether  this  temperature  dependence  is  characteristic  of  CCOs  in  general  or  is 
only  characteristic  of  our  CCD,  but  it  does  seem  to  be  reasonably  consistent 
with  the  physics  of  CCOs. 

As  we  shall  discuss  below,  the  1/f  noise  in  the  output  on-chip  FET 
amplifier  is  responsible  for  the  readout  noise.  Because  noise  in  any  of  the 
supply  voltages  couples  to  the  on-chip  amplifier,  it  is  important  to  minimize 
power  supply  noise.  We  constructed  power  supplies  using  standard  single  chip 
IC  regulators.  These  devices  have  good  noise  specifications  at  high 
frequencies  but  not  at  the  low  frequencies  at  which  we  sample  our  CCDs.  To 

remedy  this,  we  used  toroidally  wound  inductors  in  series  with  the  power  input 

to  the  regulator  chips.  The  equivalent  noise  of  our  electronics  is  5 
electrons,  which,  because  it  adds  in  quadrature  to  the  readout  noise,  only 
increases  the  readout  noise  by  a  negligible  amount. 

To  test  that  we  understand  the  origin  of  the  readout  noise,  we  have 
calculated  the  readout  noise  using  the  measured  CCD  capacitance,  the  measured 
1/f  noise  of  the  on-chip  amplifier,  and  our  sampling  function.  The  calculated 
noise  is  57  electrons  with  an  estimated  error  of  6  electrons,  and  the  measured 
readout  noise  is  64  electrons,  in  reasonable  agreement. 

The  CCD  system  is  managed  by  a  NOVA  minicomputer  which  we  use  to  download 
the  CCD  control  microprocessor  and  to  analyze  images  at  the  telescope.  The 
NOVA  is  programmed  in  FORTH  and  has  a  variety  of  routines  to  perform  simple 
picture  processing  at  the  telescope.  A  particularly  convenient  routine  that  we 
have  written  for  focussing  the  telescope  when  the  CCD  is  in  operation  displays 
a  cross-section  of  the  stellar  image  on  the  TV  monitor  and  also  displays  the 
calculated  full -width-half-maxi mum  in  pixels.  Since  the  display  is  updated 
every  second  it  greatly  simplifies  focussing. 

Image  processing  at  Michigan  is  carried  out  on  a  VAX/780  computer. 
Attached  to  the  computer  is  a  Grinnell  GMR270  which  is  capable  of  displaying 
images  on  a  color  TV  monitor  in  4096  pseudocolors. 

The  image  processing  system  that  we  use  started  out  as  a  set  of  basic 
image  processing  programs  written  at  Mt.  Stromlo  Observatory.  We  have  added 
to  the  programs  so  that  we  can  analyze  data  from  any  two-dimensional  detector 


and  produce  outputs  on  a  variety  of  output  devices.  Programs  are  now  available 
to  perform  arithmetic  functions  on  images,  to  carry  out  photometry  on  both 
stars  and  galaxies,  to  flat-field  CCD  images,  and  to  analyze  speckle 
interferometry  images. 

Stellar  photometry  can  be  easily  carried  out  on  stars  covering  a  wide 

range  of  intensities.  The  brightest  stars,  which  have  the  highest  photometric 
accuracy,  cannot  be  used  for  calibration  on  most  CCD  systems  because  for 
typical  exposure  times  they  saturate  the  system.  We  use  bright  stars  for 
calibration  by  using  short  shutter  speeds  and  by  defocussing  the  stellar  image. 
This  technique  allows  us  to  obtain  photometric  calibrations  to  better  than  one 
per  cent  on  the  best  bright  reference  stars. 

Considerable  effort  has  gone  into  writing  software  for  the  reduction  of 
two-dimensional  surface  photometry  of  elliptical  galaxies.  We  now  have 
routines  which  can  do  a  weighted  least-squares  fit  of  any  chosen 

surface-brightness  formula  to  the  surface  brightness  distribution  of  an 
elliptical  galaxy  wh^le  taking  into  account  overlapping  stars  and  galaxies, 
including  "seeing  effects"  due  to  the  earth's  atmosphere,  yielding  the  fitting 
parameters  of  the  formula  with  the  associated  errors.  Some  of  our  surface 
brightness  measurements  reach  to  beyond  30th  magnitude,  fainter  than  any 
previously  published  results. 

One  interesting  application  of  our  elliptical  galaxy  fitting  routine  was 
to  fit  a  bright  elliptical  galaxy  that  was  partially  overlapping  the  images  of 
other  faint  ellipticals  and  to  completely  subtract  out  the  image  of  the  bright 

galaxy.  This  method  could  also  be  used  to  obtain  images  of  individual  objects 

in  a  cluttered  field.  Elliptical  contours  are  a  fair  approximation  to  the 
isophotes  of  a  symetrical  object. 

A  group  of  programs  has  been  written  to  reduce  speckle  interferometric 
data.  This  data  is  obtained  from  a  64x64  area  of  pixels  closest  to  the  readout 
amplifier  to  minimize  charge  smearing  by  minimizing  the  number  of  charge 
transfers  to  the  readout  amplifier  and  to  perform  the  readout  as  quickly  as 
possible.  The  power  spectrum  of  several  hundred  images  of  a  binary  star  is 
divided  by  the  power  spectrum  of  a  corresponding  number  of  images  of  a  nearby 
calibration  star.  From  the  ratio  of  power  spectra  it  is  possible  to  obtain  the 
angular  separation,  position  angle,  and  relative  intensity  of  an  otherwise 
unresol vable  binary  system.  The  advantage  of  our  system  over  other  speckle 
interferometric  systems  is  that  we  use  no  image  intensifier,  which,  though  it 
limits  us  to  stars  brighter  than  10th  magnitude  on  a  52  inch  telescope,  takes 
advantage  of  the  high  intrinsic  stability  of  the  CCD. 

We  shall  now  discuss  scientific  research  accomplishments  during  the  period 
of  the  grant.  Below,  we  shall  discuss  scientific  progress  in  four  different 
research  areas:  global  cosmological  tests,  the  resolution  of  the  binary  star 
system  Mu-Cassiopeiae  and  a  measurement  of  the  primordial  helium  abundance,  and 
cosmological  dark  matter. 

Global  cosmological  tests  are  the  best  way  to  determine  the  future 
expansion  of  the  universe  —  Will  it  expand  forever  or  will  it  eventually 
collapse  and  fall  back  on  itself?  At  the  present  time  it  appears  possible  to 
use  CCDs  to  determine  the  answer  to  this  question.  We  have  been  working  on  two 
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cosmological  tests,  the  metric  angular  diameter  and  the  redshift-magnitude 
tests.  It  appears  to  be  possible  to  carry  out  both  tests  using  the  same  set  of 
data.  The  problems  at  the  present  time  center  on  being  able  to  reduce 
systematic  errors  in  the  collected  data  so  that  they  do  not  bias  the  results. 
At  the  present  time,  our  work  looks  promising  and  we  believe  that  we  shall  be 
able  to  separate  the  two  most  troublesome  systematics,  mergers  of  galaxies  and 
stellar  luminosity  evolution. 

We  are  in  the  process  of  writing  up  a  study  on  the  distant  cluster  of 
galaxies,  Abell  1689.  We  have  collected  photometric  data  in  two  colors  on  over 
200  elliptical  galaxies  in  the  cluster.  This  rich  cluster  of  galaxies  consists 
of  many  overlapping  images  and  has  never  been  adequately  analyzed.  Also,  many 
of  the  galactic  images  are  quite  small  and  have  been  smeared  out  significantly 
by  atmospheric  turbulence  so  that  it  is  necessary  to  correct  for  such  effects. 
We  accomplish  this  by  using  the  Images  of  stars  in  the  CCD  frame  that  surround 
the  galaxy  of  interest. 

By  making  a  plot  for  the  cluster  of  the  characteristic  size  versus 
magnitude  of  the  galaxies,  it  has  been  possible  to  find  a  range  of  galaxies 
that  show  no  dynamical  evolution.  Because  only  the  brightest  galaxies  show 
effects  of  mergers  and  only  the  faintest  galaxies  show  effects  of  tidal  . 
distortion,  there  is  a  middle  range  of  galaxies  that  reveal  none  of  these 
dynamical  problems.  However,  these  galaxies  are  affected  by  stellar  evolution. 
But,  if  only  stellar  evolution  is  present,  it  is  possible  to  deduce  the  amount 
of  it  because  the  surface  brightest  of  galaxies  should  scale  as  the  fourth 
power  of  the  redshift.  Any  deviation  from  it  indicates  stellar  evolution,  the 
amount  of  which  can  then  be  deduced.  Perhaps  the  most  Interesting  aspect  of 
this  work  from  the  point  of  view  of  other  applications  is  that  it  is  an  example 
of  treating  data  consisting  of  small,  faint  images  which  must  be  corrected  for 
atmospheric  smearing  in  a  statistically  correct  manner  to  extract  as  much 
information  as  possible. 

We  have  been  making  progress  towards  the  resolution  of  the  components  of 
the  binary  star  system  Mu-Cassiopeiae.  The  two  stars  are  separated  by  about  an 
arc  second  and  differ  by  a  factor  of  several  hundred  in  intensity  in  the 
visible.  It  is  the  close  separation  and  the  large  difference  in  intensity 
between  the  two  stars  which  has  made  the  problem  so  difficult.  With  one 
accurate  measurement  of  the  separation  between  the  components  it  is  possible  to 
determine  the  mass  of  the  stars  and  then,  with  other  information  already 
available,  it  is  possible  to  use  the  theory  of  stellar  interiors  to  determine 
the  helium  abundance  of  the  stars.  Data  has  been  obtained  on  Mu-Cas  during 
several  different  runs.  On  each  run  some  new  information  was  obtained  so  that 
the  observing  technique  and  the  method  of  analysis  gradually  improved  to  the 

extent  that  this  past  January  we  were  able  to  obtain  some  very  high  quality 

data  which  has  yielded  an  accurate  measurement  of  the  separation  between  the 
stars.  One  of  the  key  points  that  we  needed  to  learn  is  that  it  is  necessary 
to  choose  a  calibration  star  as  close  on  the  sky  to  Mu-Cas  as  possible  and  to 
also  observe  it  as  close  in  time  as  possible.  By  limiting  the  spectral 
passband  of  the  instrument  to  about  1000  Angstroms  we  found  that  it  was  not 
necessary  to  choose  a  star  that  matched  the  spectral  characteristics  of  Mu-Cas. 

A  long  term  project  on  which  we  have  been  working  has  been  to  study  the 

nature  of  the  dark  matter  in  the  halo  of  spiral  galaxies.  Evidence  to  support 


the  existence  of  dark  matter  In  halos  Is  quite  compelling.  The  rotation  curves 
of  spiral  galaxies  must  satisfy  the  equilibrium  condition  that  the  centripetal 
acceleration  necessary  for  circular  motion  be  due  to  the  gravitational 
acceleration  of  the  matter  In  the  galaxy.  However,  the  matter  contained  In  the 
visible  stars  In  galaxies  Is  not  sufficient  to  explain  the  observed 
gravitational  acceleration.  The  conclusion  Is  that  there  must  be  large  amounts 
of  unseen  mass.  The  problem  on  which  we  have  been  working  Is  to  try  to 

understand  the  nature  of  this  dark  matter. 

One  part  of  the  study  has  been  to  obtain  the  best  limits  on  the  surface 
brightness  of  the  halo  of  the  edge  on  spiral  galaxy  NGC  4565  because  it  has  a 
well  determined  rotation  curve.  Then,  with  measurements  of  the  surface 

brightness  of  the  halo,  it  is  possible  to  calculate  the  mass-to-light  ratio  of 
the  dark  matter  in  the  halo.  We  have  previously  made  observations  with  an 
Instrument  designed  to  search  for  low  surface  brightness  objects.  But,  we  have 
also  used  our  CCD  system  to  obtain  a  better  limit  on  the  surface  brightness  of 
the  halo  of  NGC  4565.  This  is  actually  a  difficult  observing  problem  for  CCDs 

because  the  halo  is  considerably  larger  than  the  field  of  view  of  the  CCD  and 

so  close  to  the  galaxy  it  is  not  possible  to  observe  the  halo  and  the 
background  sky  in  the  same  exposure.  Instead,  It  was  necessary  to  take  three 
exposures,  stepping  out  from  the  center  of  the  galaxy  with  overlap  between  the 
exposures.  With  this  technique  it  has  been  posible  to  reach  surface 
brightnesses  at  the  level  of  a  few  ten-thousandths  of  the  sky.  By  fitting  the 
observed  surface  to  the  expected  1/r  form  for  the  projected  distribution  of  the 
dark  matter,  it  is  possible  to  determine  an  upper  limit  on  the  mass-to-light 
ratio  of  the  dark  halo  matter. 


K.  Olive  (Fermi lab)  and  0.  Hegyi  have  used  this  data  to  present  a  series 
of  arguments  concluding  that  it  is  quite  unlikely  that  the  dark  halos  can  be 
composed  of  baryonic  matter.  Specifically,  a  halo  of  gas,  dust  and  rocks, 
snowballs,  low  mass  stcirs,  planets,  dead  stars  or  neutron  stars  is  shown  to  be 
unlikely.  We  briefly  summarize  the  arguments.  A  gaseous  halo  will  quickly 
collapse  and  heat  up  giving  off  more  X-rays  than  Is  consistent  with  upper 
limits  on  the  X-ray  background.  A  halo  composed  of  dust  or  rocks  would  have 
contaminated  the  gas  that  formed  stars  In  the  galactic  disk  with  an  excessively 
high  metal  abundance.  Snowballs  would  have  evaporated  very  quickly.  Low  mass 
stars  will  give  off  too  much  light  to  go  undetected  In  the  optical  observations 
discussed  above.  Planets,  since  they  may  be  expected  to  be  a  continuation  of 
the  same  mass  distribution  function  that  lead  to  stars,  are  also  excluded  when 
one  considers  the  arguments  in  detail.  Finally,  black  holes  and  neutron  stars 
are  not  allowed  because  the  ejecta  of  the  parent  star  that  leads  to  either  a 
white  dwarf  or  neutron  star  would  contaminate  with  too  high  a  metal  abundance  a 
latter  generation  of  stars  that  formed  In  the  disk. 
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The  design  and  >..!■, islrin  lun' m  a  i  li.irge  onipled  device  (CCD)  imaging  system  iiplinii/.ed  for  tlic  detec¬ 
tion  of  low  siirf;ua-  in  ighiiu  ss  ol)|ta  is  is  discussed  in  detail.  In  particular,  we  desc  ribe  the  design  of  the 
Dewar  and  llu  (cl)  elcelioiiics,  us  well  .1^  the  choice  of  CCD  operating  voltages.  cliKk  liniings,  and 
operating  tennH-iaiurc.  In  addition,  we  discuss  the  incasurenienl  ol  readout  and  spatial  noise  for  an 
RCA  SID53612  thinned  CCD  by  512  pixels).  For  this  chip,  we  find  a  readout  noise  of  64  electrons 
and  obtain  a  well  depth  in  ( mcss  of  .150  (XX)  electrons  with  essentially  no  blemishes  or  bleeding  of 
charge  on  brigtii  Mar  images  1  he  thieshold  level  for  trailing  of  charge  is  3(X)  elections  in  the  ccnier  of 
the  chip  and  is .  sseaiiallv  zeio  iit  llic  corner  of  the  chip  which  is  read  out  first  On  hioadband  exi-Hisures, 
fringing  from  niglil-sky  lines  i\  typically  seen  only  in  the  infrared  j~87(X)  A|  at  about  2%  ol  the  night - 
sky  level.  Our  twilight  Hat  field  e,\posuies  taken  on  different  nights  arc  consistent  with  each  other  at  a 
level  varying  from  one  par!  in  KXJOtX'  to  a  few  parts  in  10  000  on  the  large  scale,  and  the  zero  level  is 
stable  to  J:  10  eUetrons  over  a  uigbt  I  he  system  is  linear  over  its  entire  range,  and  we  can  currently  do 
photometry  to  better  than  0  ,  using  standard  stars. 
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I  INlKi)nLClI(i\ 

1  he  invention  of  the  eharge-eoiipleJ  ciCMce  (CCD)  by 
.Ann  iio,  Tomp.sett,  and  Smith  i  i‘)70i,  aiul  its  subsequent  de- 
V,  iopinent  over  the  past  Jecaile.  has  pimided  astrophysi- 
risi,-.  with  an  almost  ideal  photon  detector  bet  ween  the  wave¬ 
lengths  of  0.35  and  l.O  //.  A  CCD  is  .1  solid-stale  device 
'.vhich,  when  struck  by  photons,  will  produce  electron-hole 
pain;  with  high  efficieney  The  electrons  are  stored  in  a  mo- 
laic  of  electrostatic  potential  wells  formed  by  a  polysilicon 
•,  ketrode  structure  and  ion  implants,  while  the  holes  are  re¬ 
pelled  from  the  wells.  In  the  thick  version  of  the  RCA  CCD 
SID52501I,  the  photons  must  pass  through  the  polysilicon 
electrodes  whch  are  highly  .ibsorbing  in  the  CV,  whcrea.s  in 
the  thinned  CCD  (SID5.161 2l,  the  silicon  substrate  is  etched 
to  e.xtreme  thinness  so  that  the  photons  can  enter  through 
the  backside  of  the  chip,  avoiding  the  absorptive  electrode 
structure.  In  both  CCD’s,  the  entire  imaging  area  ts  light 
sensitive:  there  are  no  “dead  regions”  along  the  edges  of  pix¬ 
el-..  Also,  because  the  thick  CCD  is  many  times  more  sensi¬ 
tise  to  cosmic  rays  than  the  thinned  chip,  the  thinned  CCD 
li.as  been  preferred  for  use  as  a  photodetector  by  astrophysi¬ 
cists  when  long  integration  times  are  required. 

lo  read  out  an  im.age  stored  in  the  CCD,  the  electrode 
structure  which  helps  form  the  potential  wells  is  used  to 
transfer  the  charge  from  the  pixels  to  a  charge-sensitive  out¬ 
put  amplifier  located  on  the  chip.  In  operation,  charge  is  first 
shdtcd  vertically,  one  row  at  a  time,  to  a  horizontal  register, 
rile  charge  is  then  shifted  out  of  the  horizi  nlal  register  onc 
ec.luinn  at  a  time  to  the  on  ..  hir  amplifier  eoiisisting  of  a 
sinal!  capacitor  (  —0.5  pip  wi.ese  voltage  is  sampled  by  an 
‘.  n-chip  FET  source  followi  1  I'he  charge  111  voltage  of  the 
capacitor  is  directly  propcrii'  iial  to  the  amount  of  charge 
d.:positcd.  and,  consequently ,  i-  iiireclly  pi  .j-ni t tonal  to  the 
'I'jiVibcr  of  photons  incident  on  .t  pixel. 

The  advantages  of  a  CC  I)  are  high  qu. -11111111  elliciencv 
iSO'T  peak),  spatially  digi'i  ed  oiitpul.  geom-.l: ic  -lability, 
la 'ge  dynamic  range  i55(X!i.  and  excehent  liiieair  y  All  exist¬ 
ing  CCD's,  however,  sui'ei  froni  iionzer.'  le.idout  noise 
ranging  from  about  10  to  10  '  eleciri.ns  .and.  for  viry  low- 
light  lev-els,  a  less  than  [■'.■ib  .:  ch.irgc-Ir.m-.ler  ellieieiiey 


which  gives  rise  to  a  “threshold  effect"  which  can  compli¬ 
cate  applications  to  spectroscopy  and  astrometry.  However, 
at  higher  background  light  levels,  the  threshold  effect  is  no 
longer  a  problem,  and  at  still  higher  background  light  levels, 
the  readout  noise  is  dominated  by  shot  noise.  Blemishes  and 
fringing  (the  interference  of  light  between  the  front  and  back 
faces  of  a  CCD)  can  further  complicate  the  data  analysis. 
From  the  experiences  of  various  groups  wtho  have  built  CCD 
imaging  systems  (Low  1977  and  1981;  Dewey  and  Ricker 
1980;  Leach,  Schild.  Gursky,  Madejski,  Schwartz,  and 
Weekes  1980;  Marcus,  Nelson,  and  Lynds  1979;  Meyer  and 
Ricker  1980;  Robinson  1981;  Wright  and  Mackay  1981), 
from  our  own  experience,  and  from  inspection  of  published 
CCD  pcitures,  it  is  clear  that  the  above  problems  are  serious. 
Accordingly,  we  shall  discuss  aspects  of  the  CCD  system 
relevant  to  these  issues  in  greater  detail. 
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ii'.  1  Schematic  ilugi.m  ol'il'c  C  i.  I )  IKvvji 

11  SYSThM  t  ONIKOI 

l  iguu  I  sluiws  an  overview  ofour  CCD  system.  We  use  a 
Data  General  NOVA  3/12  niiii'eoiiiputer  to  provide  the 

■  iverall  control  of  the  CCD  sy.^lem.  I'he  selection  of  a  NOVA 
a.s  a  host  computer  may  seem  somewhat  dated  at  the  present 
time,  but  in  1978,  when  the  design  and  construction  v.as 
started  by  one  of  us  (DH),  it  was  a  reasonable  choice  and  was 
■ivailable  from  another  project.  The  NOVA  supervises  an 
.'sposure  by  loading  the  memory  of  the  CCD  control  micro- 
piocessor  located  near  the  CCD  by  opening,  timing,  and 
closing  t  he  shutter,  by  initiating  the  readout,  and  by  transfer- 
:  ing  the  data  to  the  magnetic-tape  unit  for  storage.  In  addi- 
don,  data  can  be  transferred  from  tape  to  floppy  disk  for 
processing,  to  a  TV  monitor  via  the  Lexidata  for  display,  or 
to  the  TTY.  All  the  NOVA  programming  has  been  carried 
out  in  lORTii  which  has  proven  to  be  very  convenient. 

We  sliall  now  discuss  the  sysi cm  in  more  detail.  Commun¬ 
ication  between  the  NOVA  minicomputer,  located  in  an  ob- 

■  civing  room  near  the  dome,  and  the  CCD  electronics, 
'inninted  beneath  the  Dewar  on  the  telescope,  is  via  a  serial 
hnl;  ovxr  a  50-ft  cable.  The  NOVA  transmits  (downloads) 
machine  code  over  the  cable  in  the  form  of  256  48-bit 
numbers  to  the  memory  of  an  .AM2b|0  microprixessor 
vhich  generates  the  necessary  timing  signals  for  control  of 
'  he  CCD.  We  chose  this  particular  microprocessor  because 
II  was  the  fastest  microprocessor  available  at  the  time  we 
designed  the  system.  Its  cycle  time  is  ICKJ  ns  Additional  ma¬ 
chine  code  is  used  for  control  of  the  shutter,  for  the  filter 
wheel,  and  to  adjust  voltages  on  the  CCD. 


S\MIM  D1 


Fk;  .1  View  of  the  lop  of  Ihe  opened  Dewar 


Readout  of  the  CCD,  sampling  and  digitization  of  1  he  out¬ 
put,  and  transfer  to  magnetic  tape  is  accomplished  in  5(t  s  by 
making  use  of  a  double  buffer  programmed  into  the  NOV  A's 
memory.  During  an  exposure,  the  time  remaining  in  the  inte¬ 
gration  •$  displayed  on  the  inoiiitor.  A  maximutti  of  55  full- 
frame  exposures  can  be  stored  on  a  1200-fl  reel  ){ tape. 

Through  the  T1  Silent  700  hard-copy  terminal,  a  use  can 
enter  a  variety  of  commands  into  the  system  for  displaying 
gray-scale  data,  graphing  data,  rapidly  leading  out  a 
64  X  64-pixel  region  for  focusing  Ihe  telescope,  performing 
speckle  interferometry,  doing  photometry,  and  performing 
picture  arithmetic.  (The  hard  copy  has  proven  useful  for  re¬ 
solving  discrepancies  in  the  observing  log.) 

Under  software  control,  the  u.ser  may  display  any  four  bits 
of  the  14  bits  of  digitized  intensity  information  from  any 
256x240-pixel  subsection  of  the  CCD  frame  by  loading  the 
data  into  the  Lexidata  200D  video  generator.  The  user  can 
also  graph  horizontal,  vertical,  or  diagonal  cuts  through  the 
data  to  further  examine  a  feature.  These  cuts  can  be  used  to 
calculate  the  full-width-half  maximum  (FWHM)  of  a  fea¬ 
ture,  which  is  useful  for  several  modes  of  operation,  but  most 
typically,  focusing.  In  the  focusing  mode,  a  repetitive  cycle  is 
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I.  .'.it'lishcd  in  which  a  ()4  c  i  piM  l  sc.  non  ol  ihc  (.'C'li  sui  - 
K  muling  a  star  or  a  spccir.il  line  is  liispUycil  aiui  oplinnal 
V  Ills  Ihrough  the  image  in  ihc  hoii/onl.il  .iinl  si  ilual  ilncc- 
timis  with  coiiCJitMiiuiing  I  WIIM  s  nia\  al.o  I'c  ilisplauil. 

I  lie  64  A  64-pixel  region  is  .ilso  iiscJ  tor  spci  klc  intcrlcrome- 
'•  1  > 

A  cursor  can  be  liispl.is.  il  on  the  vulco  inoinlor  to  iii.nk 
ot  I  i.l''  to  be  pholomeleicvl.  pixel  \  lines  in  he  le.nl  out,  .iinl 

■  ions  to  be  averaged.  Some  simple  piclin.-  processing  ion- 
n  .. .  are  available  so  th.it  subtle  features  in  images  mav  be 
si, idled  while  still  at  the  observatory  lo  see  if  additional  e\- 
j  •  sines  are  necessary.  Ronliiiesarcavailatile  to  perform  llal- 
l!  kiing,  dark  current  and  bias  subtract  ion,  and  to  eo  add 
images.  Fully  automatic  operation  is  possible  for  which  the 
user  Is  not  required  to  be  present;  this  is  useful  for  taking  bias 
i'..i  lies  and  dark  frames. 

iinm.  c'(  i)!)t;\V  AR 

ihc  liquid  nitrogen  Dewar  for  the  CCD  is  a  cylindrical 
.li.miber  7.75-in.  in  diaiii  and  8.5-in.  in  length  with  a  9- 
n,  diam.  faceplate  (see  Fig.  2).  It  is  designed  to  work  in  an 
upward-  or  downward-looking  position  by  inlerchaiiging 
■  Lily  the  external  fill  and  vent  lube  conneetions.  I'he  CCD 
s.ieket  is  screwed  to  a  copper  block  and  the  C'CDchip  is  held 

II  place  with  two  copper  fingers  The  copper  block  is  at- 
i.c-hed  by  means  of  stainle.ss  steel  standoffs  lo  a  1-liter  weld¬ 
ed  Stainless  steel  liquid-nil rogen  bottle.  The  nitrogen  bottle 
i.s  supported  by  a  combination  of  Kel-F  standoffs  and  by  the 
0  5  in.  diam  20-mil  stainless  steel  (ill/vent  tube.  Three  addi- 
'  nal  Kel-F  rods  connect  the  cylindrical  radiation  shield  to 
Hie  outer  Dewar  wall.  Flexure  is  undetectable  in  this  ar- 
Kuigement.  and  the  nitrogen  holding  time  is  nut  seriously 
degraded  by  our  rigid  structure;  one  tilling  lasts  for  10.5  hr. 
A  photograph  of  the  top  of  the  opened  Dewar  is  shown  in 
Fig.  3. 

IV.  THF  CCD  ELI  c  IRONICS  AND  LOGIC 

figure  4  shows  a  bltK'k  diagram  of  the  CCD  electronics, 
t  lie  digital  signal-processing  circuitry  is  enclosed  in  a  heavi¬ 
ly  shielded  and  temperature-controlled  chassis  divided  into 
separate  sections  and  connected  to  the  Dewar  by  a  triply 
shielded  1 8-in. -long  cable.  Power  supplies  are  mounted  out¬ 
side  the  box  and  power  is  brought  in  through  filtering  feed¬ 
throughs.  Although  we  initially  placed  a  preamplifier  inside 
the  Dewar,  we  found  that  greater  stability  could  be  achieved 
wuhout  it  and  at  no  incre.ise  in  readout  noise.  Our  output 
hies  level  is  presently  stable  to  ±  10  electrons  over  a  night. 

J II  the  following  p.aragraphs,  we  discuss  the  preparation 
ai.d  readout  of  the  CCD  and  trace  the  signal  through  the 
electronic  circuitry.  In  order  to  prepare  the  system  for  an 
exjKisure  and  readout,  the  NOVA  downloads  a  variety  of 
subroutines  in  machine  code  into  the  microprocessor  mem¬ 
ory  The  first  microprocessor  memory  location  is  reserved 
till  a  jump  instruction  which  moves  the  program  pointer  to 
ill.,  desired  program.  When  the  Nova  initiates  a  particular 
command,  it  only  loads  the  address  of  the  chosen  program 
mto  the  first  microprcKessor  memory  location  and  starts  the 
microprocessor.  An  exposure  and  its  preparation  consist  of 
the  execution  of  FORTH  commands  on  the  NOVA  and  mi¬ 
croprocessor  programs  which  (il  rapidly  read  out  the  CCD 
sic  times,  after  which  all  residual  charge  from  a  previous 
cx[v.isure  is  flushed  from  the  system,  liii  apply  the  correct 
voltages  for  an  integration  ol  the  CCD  electrodes,  and  final¬ 
ly,  mil  open  the  shutter  for  a  predetermined  period  of  time. 
Immediately  upon  completion  ol  an  exposure,  a  different 


jump  instruction  is  scnl  lo  the  CCD  and  it  is  eliK'ked  in  its 
normal  shiw -readout  miHle 

We  till II  now  lo  the  delails  ofshifliiig  Ihe  t  liargeoii  a  pixel 
acioss  llie  CCD  siibsl  rate  lo  the  output  amplifier  by  stepping 
voltages  on  the  hori/mital  and  vertical  electrodes,  and  lo  the 
sampling  and  digiliiration  of  the  charge  The  charge  on  a 
pixe-l  is  most  accurately  dcici  iimie'd  by  [Kiforining  a  differ 
cntial  sample  of  the  voll.igc  of  the  on-chip  CCD  cai>acilor 
before  and  after  eharge  is  iraiisici  rcd  lo  it.  In  our  system,  wc 
use  an  up-down  integrator  foi  this  task. 

Wc  shall  outline  Ihc  csscniial  steps  in  the  diffcrenlial  sam¬ 
ple,  referring  the  reader  lo  I-igs  5  -7  for  more  specific  details 
of  the  waveforms.  (Our  ditfercntial  sampler  is  based  on  a 
design  by  Marcus  cl  al.  1979f  First,  the  CCD  capacitor  is 
reset  to  a  reference  voltage,  and  that  voltage  is  sampled  by 
up-integrating,  when  switch  #2  (integrator)  is  turned  on. 
Then  the  integrator  is  turned  off  and  signal  charge  is  trans¬ 
ferred  to  Ihe  CCD  capacitor  by  driving  the  three  horizontal 
electrodes  ihrough  a  cycle.  To  perform  the  down-integra¬ 
tion,  switch  #  I  is  reversed,  and  the  integrator  is  turned  on 
again.  After  the  end  of  the  down-integration,  the  voltage  on 
the  integrator  is  digitized  lo  14  bits  and  the  integrator  reset 
with  switch  #3  in  preparation  for  the  next  cycle.  After  320 
horizontal  advances,  charge  is  shifted  by  one  row  in  the  ver¬ 
tical  direction  loading  the  horizontal  register,  and  another 
320  horizontal  shifts  are  performed. 

As  may  be  seen  in  Fig.  7,  we  use  three  distinct  voltage 
levels  on  the  verticals,  in  contrast  to  other  groups  that  use 
only  two  voltage  levels.  The  intermediate  slate  is  used  on  all 
vertical  electrodes  during  the  exposure  and  on  vertical  elec¬ 
trode  #2  during  the  horizontal  transfer.  We  found  that  the 
presence  of  an  intermediate  state  increases  the  charge-trans¬ 
fer  efficiency  and  the  depth  of  the  potential  wells. 

Including  the  timing  variables  shown  in  Fig.  8  which  show 
the  8  r's  for  the  fast  reads  used  for  reading  the  area  surround¬ 
ing  a  64  X  64-pixel  area,  there  are  a  total  of  29  timing  values. 
These  values  are  software  variables  and,  consequently,  are 
easily  changed.  There  are  nine  independent  electrode  vol¬ 
tages  which  are  set  using  hardware. 

Trimming  the  differential  sampler  makes  use  of  both 
hardware  and  software  adjustments.  The  product  of  gain 
X  time  must  be  the  same  for  both  the  up-  and  down-integra¬ 
tions  during  sampling.  Three  potentiometers  (not  shown)  in 
the  output  arms  of  the  switches  in  Fig.  5  can  be  adjusted  to 
balance  the  gain.  Final  adjustment  is  most  conveniently 
made  by  fine-tuning  the  ratio  of  the  time  spent  in  the  up  and 
down  portions  of  the  integration.  For  example,  our  time  in¬ 
terval  for  up-integration  r,  is  124,7  /zs  and  the  time  interval 
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fio.  5.  Schematic  diagram  of  the  differential  Mimpter.  In  the  depiction  of  the 
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: or  till  vn-iiitegration  r, ,  is  125.()//s.  Other  timings  arc  dis- 
e  j'Si  d  III  detail  in  Sec.  V,  on  optinii/ation.  Witli  the  system 
halaiictd,  the  noise  due  to  the  eleetronies  is  about  live  elec- 
irons  Compared  to  the  readout  noise  016,19  electrons,  the 
n  use  due  to  the  electronics  makes  a  negligible  contribution 
t  '  Ihe  overall  noise  since  these  components  add  in  quadra- 
tine 

Because  differential  sampling  is  not  widely  understotxl, 
w  e  vvish  to  emphasize  the  importance  of  keeping  the  ratio  of 
inlegratuin  time  to  transfer  time  large  when  l//noise  is  sig- 
nific.ini,  as  is  the  case  with  CCD  on-chip  amplifiers.  While 
the  contribution  to  1// noise  is  independent  of  integration 
time  if  one  neglects  the  transfer  time,  a  nonzero  transfer  time 
results  in  a  signal-to-noise  ratio  (SNR)  proportional  to 
;  1  Ai)',  where  at  is  the  ratio  of  transfer  time  to  total  integra¬ 

tion  plus  transfer  time  (Loh  1977;  Hegyi  and  Burrows  1980). 
Our  luirizontal  transfer  time  is  3.4  //s;  hence,  the  SNR  is 
degraded  by  2.7%  compared  to  an  infinitesimally  small 
transfer  time.  Ifthe  transfer  time  is  12/<sand  the  integration 
time  18  /IS  for  each  portion  of  the  differential  sample,  as  in 
some  systems,  the  degradation  is  about  44% .  It  is  interesting 
lo  cuinpare  the  SNR  obtained  w  ith  the  differential  sampling 
f  uiicti  ui  w  hich  we  use  for  our  CCD  w  ith  that  of  a  sampling 
liinciinn  optimized  for  a  combination  of  t//'and  white  noise 
iHcgyi  and  Burrows  1980;  Hegyi  and  Burrows  1984).  For 
small  signals,  an  optimal  sampling  function  could  improve 
ihc  SNR  by  —40%.  To  derive  this  result,  we  have  estimated 
that  the  white-noise  power  is  equal  to  the  l//noise  power  at 
.■bout  UK)  KHz. 


Two  problems  we  did  not  anticipate  during  the  initial  de¬ 
sign  of  the  CCD  electronics  which  deserve  mention  are  the 
generation  of  light  in  the  on-chip  amplifier  under  certain 
operating  voltages  and  the  drainage  of  charge  during  long 
exposures  from  several  rows  of  the  CCD  along  the  horizon¬ 
tal  register.  We  have  since  heard  that  the  latter  problem  has 
been  observed  by  others.  Our  solution  to  both  these  prob¬ 
lems  has  been  to  disconnect  the  on-chip  amplifier  output  and 
keep  the  phase  reset  voltage  low  during  an  exposure.  A  large 
load  resistance  for  the  CCD  is  only  a  partial  solution  to  these 
problems  because  although  it  modestly  increases  the  John¬ 
son  noise  at  the  input  of  the  first  amplifier,  it,  more  impor¬ 
tantly,  decreases  the  signal  power  received  from  the  CCD 
relative  to  the  background-noise  power.  A  lower  CCD  tem¬ 
perature  also  decreases  the  charge  leakage,  but  the  tempera¬ 
ture  is  constrained  by  factors  discussed  in  Sec.  V,  By  discon¬ 
necting  the  amplifier  output  we  are  able  to  use  a  mi^erately 
low  resistance  in  the  source  follower,  2000 12. 

V  OPTIMIZATION  OF  CCD  PERFORMANCE 

We  spent  considerable  time  and  effort  optimizing  the  per¬ 
formance  of  our  CCD.  Our  CCD  exposures  are  now  charac¬ 
terized  by  excellent  resolution,  large  dynamic  range,  and  an 
imaging  area  essentially  free  of  blemishes.  The  desired  char¬ 
acteristics  for  the  performance  of  a  CCD  system  are  sum¬ 
marized  in  Table  I.  To  achieve  this,  the  parameters  in  our 
system  listed  in  Table  II  were  designed  to  be  adjustable. 
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Fio.  8  Timing  logic  for  fast  readout. 
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Tabi  E  I  Desired  characuiisiits  of  .1 1. X  l) 

A  The  CCD 

I  l.atKe dynamic  raii|te(li>a' tc.idoiit  itotM  auil\liT|) |H)ti'nUal  vvcllsf 
1.  High  resolution 
.1.  Small  number  of  hiciiiislu-s 
4.  No  leakage  of  charge  from  v.i  II  > 
i.  l  ow  dark  current 
t).  laiw  noise  pickup 
7.  Timestabifity 
H  The  Electronics 

1 .  Lowest  pouible  noiM: 

2.  No  saturation 

.V  Perfect  differential  sampler 
4.  Time  stability 

Alter  the  system  was  constructed,  we  first  adjusted  the 
parameters  to  obtain  a  usable  picture  from  the  CCD.  Next, 
AC  began  a  long  series  of  tests  to  fine-tune  the  system.  The 
jiiocess  of  optimization  is  complicated  by  the  fact  that  a  sin¬ 
gle  control  parameter  can  affect  more  than  one  CCD  charac- 
leiistic  and,  consequently,  adjusting  one  control  parameter 
may  necessitate  readjusting  another.  These  interactions, 
however,  are  not  serious,  and  in  most  cases  the  CCD's  char- 
ac  let  istics  are  unchanged  by  a  small  variation  in  the  value  of 
a  parameter. 

.A  large  dynamic  range  is  achieved  by  optimizing  the  CCD 
so  that  the  potential  wells  can  be  filled  to  their  maximum 
capacity.  With  a  readout  noise  of  about  64  electrons  for  the 
RCA  thinned  chip,  a  dynamic  range  of  about  SSOO  can  be 
atiuined.  Whereas  the  readout  noise  is  mostly  determined  by 
the  I// noise  spectrum  of  the  on-chip  amplifier  and  is  sensi¬ 
tive  to  timing  ratios  in  the  differential  sampling  function,  the 
ability  of  the  potential  wells  to  hold  charge  is  strongly  deter¬ 
mined  by  the  voltages  on  the  horizontal  and  vertical  transfer 
electrodes. 

In  order  to  determine  a  first  approximation  to  the  operat¬ 
ing  potentials  of  the  electrodes  and  so  obtain  a  usable  pic¬ 
ture,  we  exposed  the  CCD  to  light  of  an  intensity  and  dura¬ 
tion  such  that  the  wells  would  be  filled  to  their  maximum 
level.  The  empirically  determined  voltage  ranges  to  prevent 
bleeding  of  charge  are  given  inTable  III.  The  voltage  limits 
were  less  stringent  at  lower  exposure  levels;  hence,  adjust¬ 
ment  at  a  high  exposure  level  guarantees  a  large  dynamic 
range. 

Having  determined  a  range  of  electrode  potentials  which 
gave  the  largest  dynamic  range,  the  next  step  was  to  deter¬ 
mine  the  potentials  whch  gave  the  best  resolution.  Because 
the  resolution  is  adversely  affected  by  the  number  of  trans¬ 
fers  of  a  charge  packet  before  it  is  deposited  onto  the  on-chip 
capacitor  of  the  CCD  and  by  low  light  levels  (the  threshold 
effect),  we  achieved  the  greatest  sensitivity  for  improvement 
of  resolution  by  imaging  a  faint  point-like  source  (a  pinhole 
in  aluminum  foil  whose  image  was  larger  than  one  pixel) 
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onto  the  last  area  of  the  CCD  to  be  read  out.  The  resolution 
was  most  sensitive  to  the  potential  on  horizontal  electrode 
#3.  By  making  the  voltage  swing  during  transfer  larger  on 
this  electrode  than  on  those  of  the  other  horizontal  elec¬ 
trodes,  i.e.,  —  7  to  4-  6  V  compared  to  —  5  to  -t-  5  V,  the 
resolution  was  noticeably  improved.  The  potential  on  hori¬ 
zontal  electrode  #2  had  less  of  an  effect,  but  when  the  vol¬ 
tage  swing  was  set  to  only  --  4.5  to  -b  4.5  V,  defects  ap¬ 
peared  on  the  screen  and  the  resolution  deteriorated. 

The  resolution  is  also  strongly  affected  by  the  time  spent 
transferring  charge.  In  general,  a  longer  transfer  time  gives 
higher  resolutioo,  but  at  the  expense  of  increased  blemishes. 
In  addition,  the  vertical  timing  values  and  the  vertical  elec¬ 
trode  potentials  are  weakly  coupled.  Referring  to  the  timing 
logic  diagram  in  Fig.  6,  we  found  that  the  horizontal  transfer 
times  r„  and  r,  had  little  or  no  effect  upon  either  the 
resolution  or  the  blemishes.  Making  T-,  longer  (up  to  1.2  fis) 
improved  the  resolution;  however,  pulse  feedthrough  on  the 
chip  increased  the  readout  noise.  A  good  compromise  for  Ty 
was  0.5  /IS.  Making  r«  longer  also  improved  the  resolution, 
but  at  the  expense  of  increased  blemish.  Its  value  was  set  to 
0.3  fis. 

For  the  vertical  transfer  times  (see  Pig.  7),  increasing  t„, 
and  r,g  improved  the  resolution  slightly,  but  also  increased 
the  number  of  vertically  oriented  blemishes.  Increasing  r,7 
and  r,9  improved  the  resolution,  but  also  increased  the  num¬ 
ber  of  blemishes,  though  Tyu  had  no  effect  on  resolution  and 
blemishes.  By  making  the  high  state  voltage  on  vertical  elec¬ 
trode  #3  more  negative,  the  number  of  blemishes  could  be 
reduced,  but  at  the  expense  of  reduced  well  capacity  (ob¬ 
served  as  bleeding  of  the  image  at  high  exposure  levels).  The 
voltage  levels  at  which  the  number  of  blemishes  were  re¬ 
duced  and  at  which  bleeding  appeared  were  found  to  be  af¬ 
fected  by  the  vertical  transfer  times.  The  best  compromise 
for  electrode  potentials  is  listed  in  Table  IV.  The  effect  of 
temperature  on  resolution  is  discussed  at  the  end  of  this  sec¬ 
tion. 

When  the  CCD  is  operated  in  its  focusing  mode,  a  64 x  64- 
pixel  area  is  read  out  with  the  same  electrode  voltages  and 
timing  values  as  for  the  full  frame,  but  the  surrounding  pixels 
are  read  out  at  a  much  faster  rate.  The  time  intervals  rji 
through  Tyt  control  this  fast  readout  and  the  timing  logic  is 


Table  III.  Electrode  potential]  for  which  bleeding  of  charge  did  not  occur. 


Electrode  number 

High  state 

Loaraiaw 

IntermediMeiMe 

Horizontal  1 

3.5<I'<6  5 

-5.9<»'<  -1.7 

Horizontal  2 

-3.0<I' 

y<  -  3.0 

Horizontal  3 

l.7<  I' 

-7.9<T<  -3.4 

Vertical  1 

-0.8<  E 

y<  -  g.o 

y<  -  2.4 

Vertical  2 

y<  -  2  0 

y<  -  8.0 

-4.0<»'<a0 

Vcitical  3 

-o.6<  y 

y<  -8.0 

y<  - 1.0 

Note  to  Table  III 


Only  the  vertical  electradei  have  an  ialcnncdiMe  Mate. 
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shown  in  Fig.  8.  We  found  that  the  readout  noise  is  slightly 
dependent  on  r22  and  rjg  and  that  values  of  4.2/rs  give  about 
a  reduction  in  noise  compared  to  0.3  /rs. 

We  also  operate  the  CCD  in  a  mode  called  “speckle” 
which  is  designed  for  short  exposures  of  the  highest  possible 
resolution.  This  mode  is  similar  to  the  focus  mode  described 
above,  except  that  the  64  X  64-pixel  area  is  placed  in  the  up¬ 
per  left  of  the  CCD  frame  (the  Hrst  area  to  be  read  out),  and 
hence  undergoes  the  smallest  number  of  transfers  before  be¬ 
ing  sampled.  In  this  mode,  32  consecutive  frames  are  record¬ 
ed  on  tape.  We  were  able  to  increase  r,  through  so  that  the 
total  horizontal  transfer  time  was  6  /is  and  to  increase 
through  T20  to  1.2/is  without  producing  blemishes.  The  im¬ 
provement  in  resolution  was  11%  as  determined  by  measur¬ 
ing  the  FWHM  of  a  point-like  light  source  imaged  onto  the 
CCD.  Fur  our  full-frame  format  and  the  long  exposures  re¬ 
quired  by  astronomical  objects,  the  resolution  is  as  good  as  in 
the  speckle  mode;  this  is  because  the  longer  exposure  times 
allow  a  sky  background  to  build  up  which  eliminates  the 
threshold  eflfect. 

The  threshold  effect,  which  is  eliminated  with  a  back¬ 
ground  of  about  300  electrons,  was  investigated  by  imaging 
faint  point-like  images  on  an  adjustable  background.  In  the 
upper  left  comer,  however,  there  is  no  measurable  threshold 
effect  The  fast  read  r’s,  and  Tjs  were,  in  addition, 

increased  to  1.2,  0.7,  and  0.6 /is  in  order  to  improve  the 
uniformity  of  the  bias  exposure  (zero  incident  light,  zero  ex¬ 
posure  time). 

Several  schemes  for  the  timing  and  switching  logic  were 
tried.  Our  first  arrangement  (based  on  the  RCA  data  sheets) 
used  four  fewer  time  intervals  for  the  horizontal  transfer.  A 
shorter  total  time  for  transfer  was  expected  to  give  lower 
readout  noise,  as  explained  in  Sec.  IV.  With  this  switching 
logic,  the  electrodes  have  different  potentials  before  and 
after  charge  transfer  and,  because  of  feedthrough  effects  in 
the  CCD,  an  offset  voltage  developed  which  caused  satura¬ 
tion  in  the  amplifier  chain.  Because  the  time  between  inte- 
pations  is  constrained  by  its  effect  on  resolution  and  blem¬ 
ishes.  the  \/f  noise  and  the  white  noise  can  be  made  smaller 
by  lengthening  the  integration  times  (the  integration  times, 
r ,  and  r , , ,  also  control  the  gain  of  the  amplifier  chain).  Other 
optimizations  of  the  timing  included  determining  the  best 
values  to  avoid  feedthrough  and  pickup  of  timing  signals  by 
the  integrator.  In  particular,  r,,  and  r,„  should  be  slightly 
greater  than  1 .2/<s.  The  best  values  for  the  timings  are  listed 
in  Table  V. 

We  found  that  the  chip  temperature,  in  addition  to  affect¬ 
ing  the  loss  of  charge  near  the  horizontal  register  mentioned 
in  Sec.  IV,  also  controls  the  quantum  efficiency  and  the  reso¬ 
lution.  For  each  degree  drop  in  temperature,  the  quantum 
efficiency  decreases  by  about  0.3%  (tested  at  a  wavelength  of 
5500  A  and  between  107  and  173  K).  The  resolution  was 
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found  tO’  be  best  at  IS2  K,  with  the  FWHM  of  a  point-like 
image  being  41%  larger  at  10  deg  higher  or  lower  than  the 
optimal  temperature  (see  Fig.  9).  Note  that  the  minimum 
I%HM  in  Fig.  9  is  constrained  by  the  size  of  the  point-like 
image  we  used.  Although  we  can  obtain  a  FWHM  of  unity, 
e.g.,  on  spearal  emission  lines,  the  actual  value  depends  on 
the  placement  of  the  image  relative  to  a  pixel’s  boundary. 
Here  we  avoided  that  problem  by  using  a  larger  image.  We 
estimate  the  accuracy  of  our  temperature  measurements  to 
be  approximately  thm  d^rees. 

Prwise  control  of  the  CX^D  temperature  is  important  not 
only  for  assuring  that  one  is  always  within  the  optimum 
range  of  operation,  but  also  to  stabilize  the  bias  level  at  the 
input  of  the  amplifier  chain.  For  each  degree  change  in  tem¬ 
perature  of  the  CCD,  the  input  bias  (equal  to  the  voltage 
drop  across  the  CCD  load  resistor)  chaises  by  3. 1  mV.  As 
we  mentioned  in  Sec.  IV,  the  output  bias  level  is  made  rela¬ 
tively  insensitive  to  this  change  by  carefully  balancing  the 
differential  sampler  by  adjusting  r,.  We  stabilize  the  tem- 
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Fig.  9.  Dependence  of  resolution  on  CCD  tem- 
penture. 
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[xiuture  of  the  CCD  by  a  viiy  liigli-guin  fcLilback  Kk)p 
w  liif  li  drives  a  heating  element,  two  resistors  iiiountcd  inside 
a  eopper  blcK'k  onto  which  the  CCD  is  attuehed.  Operation 
at  an  average  power  input  of  about  I  W  gives  some  addi¬ 
tional  sensitivity  (depending  on  the  square  root  of  the  iwwer) 
and  does  not  seriously  decrease  the  holding  time  lor  liquid 
nitrogen  in  our  Dewar. 


VI  DETERMINAIION  Ot  (X  I)  NOISE 

1  he  readout  and  spatial  noise  of  a  CCT)  arc  of  importance 
because  the  former  determines  the  niiniinum  delectable  sig¬ 
nal  for  situations  with  a  small  miniber  of  collected  photons 
|x-r  pixel  such  as  spectroscopy  and  speckle  interferometry, 
w  hile  the  latter  sets  the  limit  for  a  large  number  of  collected 
photons  such  as  the  detection  of  faint  features  superimposed 
on  a  night-sky  background.  Although  spatial  noise  can  be 
corrected  by  dividing  by  a  flat-tield  exposure,  the  variation 
in  quantum  efficiency  with  wavelength  from  pixel  to  pixel  is 
sufficiently  large  that  differences  between  the  spectrum  of 
the  signal  and  that  of  the  flat  field  become  important  at 
roughly  the  level  of  a  few  parts  in  1000,  for  wavelengths  near 
bright  night-sky  lines. 

The  measurement  of  readout  noi.se  can  be  done  in  two 
ways.  1  he  methods  start  from  the  basic  equation 

+eri„  (1) 

where  a,,  is  the  total  noise,  a,,,  the  readout  noise,  a,  the 
shut  noise,  and  a,  ,  the  spatial  noise.  Alt  quantities  are  on  a 
per  pixel  basis  and  are  measured  in  units  of  electrons,  if  we 
define  the  gain  g  as  ADU/n,,  where  ADU  represents  ob¬ 
served  analog-to<iigital  units  and  n,  is  the  corresponding 
number  of  electrons,  and  define  F,,  the  fractional  spatial 
noise,  as  a,  ,/n,,  then  because  a,  =  {n,  we  find  the  total 
noise  in  ADU  to  be 

‘^.XOU  =  <^.AD0  +  ff-ADU  -E  f  JADU^  (2) 

The  first  method  makes  use  of  a  least-squares  fit  of  the 
observed  of, ado  versus  ADU  data  to  the  above  second-order 
equation.  It  is  not  difficult  to  show  that  each  point  must  be 
weighted  by 

K  -  ,3, 

2<0f.AD0>^’ 

where  is  the  number  of  pixels  in  each  exposure  and  is 
the  number  of  repeated  ob^rvations  of  each  point.  If  infor¬ 
mation  about  the  spatial  noise  is  not  needed,  some  computa¬ 
tion  can  be  saved  by  fitting  to  a  first-order  equation  after 
removing  the  spatial  noise  from  the  data.  There  is  an  advan¬ 
tage  to  this  second  method  in  that  any  nonlinearities  become 
readily  apparent.  The  spatial  noise  is  removed  by  averaging 
many  exposures  and  computing  the  mean  error  about  the 
mean  value  in  ADU  for  each  pixel.  A  variation  of  the  second 
method  and  our  chosen  approach  removes  the  spatial  noise 
by  taking  only  two  exposures.  It  can  be  shown  that  the  best 
estimator  of  the  total  noise  minus  the  spatial  noise,  of ,  adu  > 


2;  (ADU,  ADU-lf 

'A  ..aOU  =  - Z - ’ 

where  the  sum  is  over  all  pixels  in  the  area.  We  expose 
twice  on  an  area  of  SOX  SO  pixels.  A  sample  plot  from  which 
the  gain  and  readout  noise  can  be  obtained  is  given  in  Fig.  1 0. 
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Flo.  to.  DetermineUon  of  CCD  gain  and  readout  noise,  o,  ,  aku  the 
total  noise  less  the  spatial  noise  measured  in  ADU. 


In  this  case,  the  reciprocal  of  the  gain  1/g  is  21.3  electrons/ 
ADU  and  the  readout  noise  is  63.9  electrons. 

As  a  check  to  verify  that  we  understand  the  origin  of  the 
readout  noise,  we  can  calculate  the  noise  using  the  measured 
CCD  capacitance  and  the  measured  \/f  noise  of  the  on-chip 
amplifier.  By  measuring  the  current  into  the  CCD  through 
RD,  the  voltage  across  the  load  resistor  as  the  CCD  u  beiiig 
read  out,  and  our  cycle  time,  we  obtain  a  value  for  the  on- 
chip  c^Mcitance  of  0.62  pF.  To  determine  the  noise  spec¬ 
trum  we  measured  the  noise  of  the  (X^  output  mospet 
from  200  to  2S  000  Hz  with  a  calibrated  noise-spectnun  ana¬ 
lyzer.  Over  this  range,  the  noise  power  spectrum  closely  fol¬ 
lows  a  \/f  law  of  the  form  F^/v,  where  F„  ia  230(^Vp.  The 
measured  exponent  of  vis  1.12  ±0.05.  Estimating,  as  in  Sec. 
IV,  that  the  white-noise  power  equals  the  X/fndac  power  at 
100  kHz,  we  calculate  (Loh  1977;  Hegyi  and  Burrows  1980) 
the  readout  noise  to  be  36.7  ±  6  electrons.  The  error  repre¬ 
sents  an  estimate  of  our  overall  measuring  accuracy.  The 
measured  and  cafcubtcd  readout  noite  are  in  reasooaUe 
agreement. 

By  measuring  the  rendout  noise  as  a  flinction  of  pixd  loca¬ 
tion  on  the  CCD,  we  find  that  the  readout  noise  is  position 
dependent.  Figure  lOsbowsdatatakenintheupper-lefloor- 
ner  of  the  chip,  but  in  other  locations  the  shape  ^the  curve  is 
no  longer  straight.  For  these  areas,  the  measured  total  none 
at  intermediate  level  exposures  (up  to  200  000  electrons)  is  a 
few  percent  higher  than  the  shot  noise.  Hus  is  most  likdy 
due  to  the  greater  contribution  of  transfer  noise  and  tnqiping 
noise,  both  of  which  depend  on  the  square  root  of  the  total 
number  of  transfers  (Beynon  and  Lirab  1980).  At  higher 
leveb  (up  to  350  000  electrons,  a  limit  set  by  our  ADC),  we 
observe  the  same  noise  at  all  locations.  Also,  when  no  light  is 
incident  on  the  CCD,  the  noise  is  essentially  the  same  every¬ 
where. 

A  measure  of  the  spatial  noise  as  a  function  of  wavelength 
and  sample  size  was  obtained  by  computing  the  standard 
deviation  of  a  square  area  of  a  twilight-sky  flat-field  exposure 
as  a  function  eff  the  length  of  an  edge  of  the  square  after 
removing  the  contribution  due  to  readout  and  ^ot  noise. 
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I'  i(j.  1 1 .  CCD  ipMial  IKMM  as  a  futicllon  (if  wavelength  and 
sample  size. 


From  the  plot  of  this  data  shown  in  Fig  II,  we  see  that  the 
L'uriections  for  spatial  noise  become  increasingly  important 
at  longer  wavelengths  and  larger  pixel  sample  sizes.  Flat- 
field  exposures  taken  of  the  dawn  and  twilight  skies  on  dif¬ 
ferent  nights  are  consistent  with  each  other  on  the  large  scale 
(the  upper  third  compared  with  the  lower  third  of  the  frame) 
from  one  part  in  100000  to  a  few  parts  in  10000  peak  to 
peak. 

A  complete  correction  of  an  exposure  for  spatial  noise  and 
zero  level  is  carried  out  on  a  pixel-by-pixel  basis  using  the 
expression 


[(exposure  —  biasg)  —  (dark  frame  ■ 
(flat-field  —  bias  g- 1 


■  biaspllscale 


where  the  dark  frames  and  flat  fields  are  averages  of  many 
exposures  and  the  scale  is  a  normalization  constant  which 
could  be  set  equal  to  the  average  intensity  of  a  pixel  for  the 
average  of  the  flat-field  exposures.  Each  processed  exposure 
must  have  a  constant  bias  value  for  the  exposure,  dark  frame, 
and  flat  field  (biaSs,  biasg,  and  bias^ ,  respectively)  subtract¬ 
ed;  in  our  system,  these  are  computed  from  an  average  of  a 
30x  SO-pixel  area  which  is  read  out  immediately  following 
readout  of  the  full  frame.  Dark  frames  for  exposure  times 
from  zero  to  more  than  43  min  show  structure  on  the  large 
scale  of  about  40  electrons  peak  to  peak.  Because  this  can  be 
as  much  as  1%  of  the  sky  level,  it  is  important  to  subtract  this 
structure  from  an  exposure.  At  our  operating  temperature, 
the  dark  current  contributes  about  25  electrons  per  hour  per 
pixel  on  average,  with  small  pixel-to-pixel  variations. 

CCD's  are  quite  sensitive  to  background  low-frequency 
electromagnetic  radiation.  We  found  that  a  small  ripple  vol¬ 
tage  on  the  shutter  power  supply  filled  pixels  with  charge  at  a 
rate  several  hundr^  times  larger  than  the  normal  dark  cur¬ 
rent.  After  eliminating  the  ripple  so  that  the  shutter  voltage 
was  purely  dc,  even  under  load,  this  problem  disappeared. 
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VII  PIcniRi:  I’RIK  KSSINti 

Our  picture  processing  is  done  on  a  VAX/7W)  conipuier 
running  under  VMS.  Images  me  displaynl  on  a  (irmnell 
(iMR270.  TheOrinnell  with  joysiick-coiilrollcd  quad  cur- 
Mtrs  and  zoom  and  pan  enable  us  lo  display  images  in  4tW6 
pseudtKolors.  Two  VTIOO  lenninals  with  graphics,  a  I'rin- 
ironix.  a  rcktronix  4M2  flai  tx-d  plollcr,  .-tnd  the  Cirinnel  arc 
each  capable  of  displaying  plots. 

We  have  written  most  of  our  present  software  for  the  re¬ 
duction  of  two-dimensional  surface  photometry  of  galaxies. 
The  most  notable  subset  of  these  routines  can  do  a  weighted 
least-squares  fit  of  any  chosen  surface-brightness  formula  to 
an  elliplical  galaxy's  light  distribution,  while  taking  into  ac¬ 
count  overlapping  stars  and  galaxies,  including  seeing  ef¬ 
fects,  yielding  the  fitting  parameters  of  the  formula  with 
their  errors.  We  are  currently  writing  programs  to  handle 
two-dimensional  spectra.  Some  programs  which  were  writ¬ 
ten  at  Mount  Stromlo  Observatory,  Australia  by  Greg 
Quinn,  Reet  Valek,  aiul  Dennis  Wame  do  fast  data  transfer 
to  and  from  disk  in  a  format  called  Standard  Astronomical 
Data  (SAD).  A  320  X  3 1 2  array  of  real*4  data  can  be  trans- 
fered  in  about  3  s.  Copies  of  the  Mount  Stromlo  programs 
were  transmitted  to  us  by  Susan  Simpkin.  We  have  given  the 
acronym  MIIPS  (Michigan  Image  Processiag  System)  to  our 
picture-processing  software. 

An  example  of  the  quality  of  our  CCD  images  is  shown  in 
Fig.  12.  It  is  a  30-inin  exposure  of  the  field  surroimdiagNGC 
4874  in  the  Coma  cluster  taken  through  a  1000-A  bandpass 
filter  centered  at  3000  A.  This  exposure,  for  which  the  entire 
frame  is  shown,  has  not  been  ^t-fielded  or  edited  in  any 
way,  yet  there  is  no  evidence  of  fringing  or  of  charge  tradiag 
due  to  imperfect  charge-transfer  efficiency.  Tlicre  is  one  mi¬ 
nor  blemish  in  two  ai^iaccnt  columns  estcmfeig  over  34 
rows.  Although  the  bright  star  has  saturated  the  ADC  (more 
than  about  330  000  eiections/pixel)  in  43  pixds,  the  star  im¬ 
age  appears  unbroadened  with  no  bleeding  of  charge,  and 
even  diffraction  spikes  are  visible. 

In  other  exposures,  we  have  found  fringing  at  less  than  2% 
of  the  sky  level  in  the  I  band,  centered  at  8700  A,  but  at  3300 
A,  only  one  exposure  out  of  30  has  shown  fringing.  It  was 
observed  to  be  at  a  level  of  less  than  1%  of  the  sky  level. 

We  have  tested  the  linearity  of  the  system  in  three  ways. 
By  exposing  on  a  constant  uniform  source  of  light  and  then 
plotting  the  observed  signal  in  ADU  against  exposure  time, 
we  see  no  departure  from  linearity  from  our  shortest  expo¬ 
sures  (less  than  10  ms)  to  the  highest  level  our  ADC  can 
record(16  383  ADU or  348  938  electrons).  By  measuring  the 
total  light  in  the  images  of  calibration  stars  observed  at  a 
variety  of  exposure  times  and  for  a  variety  of  magnitudes,  we 
obtain  a  caUlnation  curve  with  a  standard  deviation  for  a 
single  observation  of  0^0073  (0.73%).  It  is  interesting  to  note 
that  we  can  calibrate  on  stars  as  bright  as  Vega.  The  exposure 
times  must  be  kept  below  1  s  and  the  images  defocus^. 

A  more  stringent  test  of  linearity  is  defined  on  the  scale  of 
a  (uxel.  That  is,  if  at  high  levels  of  illumination,  charge  drains 
to  neighboring  pixels  but  the  total  charge  is  conserved,  the 
above  two  tests  will  indicate  linearity  although  it  is  violated 
on  a  per  pixel  basis.  By  measuring  the  FWHM  of  point-like 
images  as  the  illumination  was  increased  from  very  low  lev¬ 
els  up  to  saturation  of  our  ADC,  we  observed  no  increase  in 
the  FWHM  of  the  image. 

As  an  example  of  image  processing  with  our  system.  Pig. 
13  shows  the  |»ocessed  field  of  Fig.  12  with  the  bd^ueet 
galaxy  modeled  by  our  programs  and  subtracted  from  the 
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Fig.  13  are  about  23‘:‘6  ibased  on  ms,**,  =;  0.0  for  Vega|.  A 
similar  modeling  and  subtiaclion  for  an  exposure  of  M87 
reveals  globular  clusters  to  the  very  center  of  the  galaxy. 
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